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(57) Abstract: Methods of predicting fracture 
risk of a patient include: obtaining an image 
of a bone of the patient; determining one or 
more bone structure parameters; predicting a 
fracture line with the bone structure parameter; 
predicting a fracture load at which a fracture 
will happen; estimating body habitus of the 
patient; calculating a peak impact force on the 
bone when the patient falls; and predicting a 
fracture risk by calculating the ratio between 
the peak impact force and the fracture load. 
Inventive methods also includes determining 
the effect of a candidate agent on any subject's 
risk of fracture. 
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SYSTEM AND METHOD OF PREDICTING FUTURE FRACTURES 



CROSS-REFERENCE TO RELATED APPLICATIONS 

[0001] This application claims the benefit of Provisional Application No. 
60/610,447, filed September 16, 2004. The present application incorporates the 
disclosure of this provisional application by reference. 

BACKGROUND OF THE INVENTION 
FIELD OF THE INVENTION 

[0002] This invention relates to using imaging methods for predicting fracture risk 
and/or location based on radiographs. 

DESCRIPTION OF THE RELATED ART 

[0003] Osteoporosis is among the most common conditions to affect the 
musculoskeletal system, as well as a frequent cause of locomotor pain and disability. 
Osteoporosis can occur in both human and animal subjects (e.g. horses). 
Osteoporosis (OP) occurs in a substantial portion of the human population over the 
age of fifty. The National Osteoporosis Foundation estimates that as many as 44 
million Americans are affected by osteoporosis and low bone mass. In 1997 the 
estimated cost for osteoporosis related fractures was $13 billion. That figure 
increased to $17 billion in 2002 and is projected to increase to $210-240 billion by 
2040. Currently it is expected that one in two women over the age of 50 will suffer an 
osteoporosis-related fracture. 

[0004] Imaging techniques are important diagnostic tools, particularly for bone 
related conditions such as osteoporosis. Currently available techniques for the 
noninvasive assessment of the skeleton for the diagnosis of osteoporosis or the 
evaluation of an increased risk of fracture include dual x-ray absorptiometry (DXA) 
(Eastell et al. (1998) New EnglJ. Med 338:736-746); quantitative computed 
tomography (QCT) (Carm (1988) Radiology 166:509-522); peripheral DXA (pDXA) 
(Patel et al. (1999) J Clin Densitom 2:397-401); peripheral QCT (pQCT) (Gluer et al. 
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(1997) Semin Nucl Med 27:229-247); x-ray image absorptiometry (RA) (Gluer et al. 
(1997) Semin Nucl Med 27:229-247; and U.S. Patent No. 6,246,745); and 
quantitative ultrasound (QUS) (Njeh et al. "Quantitative Ultrasound: Assessment of 
Osteoporosis and Bone Status", 1999, Martin-Dunitz, London England; WO 
9945845; WO 99/08597; and U.S. Patent No. 6,077,224 which is incorporated herein 
by reference in its entirety). 

[0005] DXA of the spine and hip has established itself as the most widely used 
method of measuring bone mineral density (BMD). Tothill, P. and D.W. Pye, (1992) 
Br J Radiol 65:807-813. The fundamental principle behind DXA is the measurement 
of the transmission through the body of x-rays of 2 different photon energy levels. 
Because of the dependence of the attenuation coefficient on the atomic number and 
photon energy, measurement of the transmission factors at 2 energy levels enables the 
area densities (i.e., the mass per unit projected area) of 2 different types of tissue to be 
inferred. In DXA scans, these are taken to be bone mineral (hydroxyapatite) and soft 
tissue, respectively. However, it is widely recognized that the accuracy of DXA scans 
is limited by the variable composition of soft tissue. Because of its higher hydrogen 
content, the attenuation coefficient of fat is different from that of lean tissue. 
Differences in the soft tissue composition in the path of the x-ray beam through bone 
compared with the adjacent soft tissue reference area cause errors in the BMD 
measurements, according to the results of several studies. Tothill, P. and D.W. Pye, 
(1992) Br J Radiol, 65:807-813; Svendsen, O.L., et al., (1995) J Bone Min Res 
10:868-873. Moreover, DXA systems are large and expensive, ranging in price 
between $75,000 and $150,000. 

[0006] Quantitative computed tomography (QCT) is usually applied to measure 
the trabecular bone in the vertebral bodies. Cann (1 988) Radiology 1 66:509-522. 
QCT studies are generally performed using a single kV setting (single-energy QCT), 
when the principal source of error is the variable composition of the bone marrow. 
However, a dual-kV scan (dual-energy QCT) is also possible. This reduces the 
accuracy errors but at the price of poorer precision and higher radiation dose. Like 
DXA, however, QCT are very expensive and the use of such equipment is currently 
limited to few research centers. 
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[0007] Quantitative ultrasound (QUS) is a technique for measuring the peripheral 
skeleton. Njeh et aL (1997) Osteoporosis Int 7:7-22; and Njeh et aL, Quantitative 
Ultrasound: Assessment of Osteoporosis and Bone Status, 1999, Martin Dunitz, 
London, England.. There is a wide variety of equipment available, with most devices 
using the heel as the measurement site. A sonographic pulse passing through bone is 
strongly attenuated as the signal is scattered and absorbed by trabecular Attenuation 
increases linearly with frequency, and the slope of the relationship is referred to as 
broadband ultrasonic attenuation (BUA; units: dB/MHz). BUA is reduced in patients 
with osteoporosis because there are fewer trabeculae in the calcaneus to attenuate the 
signal. In addition to BUA, most QUS systems also measure the speed of sound 
(SOS) in the heel by dividing the distance between the sonographic transducers by the 
propagation time (units: m/s). SOS values are reduced in patients with osteoporosis 
because with the loss of mineralized bone, the elastic modulus of the bone is 
decreased. There remain, however, several limitations to QUS measurements. The 
success of QUS in predicting fracture risk in younger patients remains uncertain. 
Another difficulty with QUS measurements is that they are not readily encompassed 
within the WHO definitions of osteoporosis and osteopenia. Moreover, no 
intervention thresholds have been developed. Thus, measurements cannot be used for 
therapeutic decision-making. 

[0008] There are also several technical limitations to QUS. Many devices use a 
foot support that positions the patient's heel between fixed transducers. Thus, the 
measurement site is not readily adapted to different sizes and shapes of the calcaneus, 
and the exact anatomic site of the measurement varies from patient to patient. It is 
generally agreed that the relatively poor precision of QUS measurements makes most 
devices unsuitable for monitoring patients' response to treatment. Gluer (1997) J 
Bone Min Res 12:1280-1288. 

[0009] Radiographic absorptiometry (RA) is a technique that was developed many 
years ago for assessing bone density in the hand, but the technique has recently 
attracted renewed interest. Gluer et al. (1997) Semin Nucl Med 27:229-247. With this 
technique, BMD is measured in the phalanges. 
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[001 0] Furthermore, current methods and devices do not generally take into 
account bone structure analyses. See, e.g., Ruttimann et al. (1992) Oral Surg Oral 
Med Oral Pathol 74:98-1 10; Southard & Southard (1992) Oral Surg Oral Med Oral 
Pathol 73:751-9; White & Rudolph, (1999) Oral Surg Oral Med Oral Pathol Oral 
Radiol Endod 88:628-35. 

[001 1] BMD does not accurately predict the presence of osteoporotic fracture. See, 
e.g., Riggs et al (1982) J Clin Invest 70:716-723; Krolner, B. and S.P. Nielsen (1982) 
Clin Set 62:329-336; Ott et aL (1987) J Bone Miner Res, 2:201-210; and Pacifici et 
al. (1987) J Clin Endocrinol Metab, 64:209-214. While BMD is correlated with long- 
term fracture risk in population based studies (Kains (1994) Osteoporosis Int 4:368- 
38 1), it cannot take into account factors that vary from patient to patient and that are 
major determinants of individual failure load and resultant fracture (Hayes, W.C. and 
M.L. Bouxsein, Biomechanics of cortical and trabecular bone: Implications for < 
assessment of fracture risk, in Basic Orthopaedic Biomechanics, V.C. Mow and W.C. 
Hayes, Editors, 1997, Lippincott-Raven Publishers: Philadelphia, p. 69-11 1; 
Kroonenberg et al. (1995) J Biomech Eng. 1 17(3):309-318; Kroonenberg et al. (1996) 
Biomechanics 29(6): 807-8 11; and Robinovitch et al. (1991) J Biomech Eng. 113:366- 
374). These factors include bone architecture and structure, bone morphology, and 
biomechanical loading and impact load. Indeed, patients receiving osteoclast 
inhibiting, anti-resorptive drugs show remarkable reductions in incident osteoporotic 
fractures by 60-65% but only small changes in BMD on the order of 4.0-4.5% 
(Reginster et al. (2000) Osteoporosis Int, 1 1(1):83-91; and Harris et al. (1999) JAMA 
14:1344_1352) 3 strongly indicating a significant discrepancy between clinical 
outcomes and BMD measurements of bone health. 

[001 2] Thus, there remains a need for compositions and methods for predicting 
fracture risk. 

SUMMARY OF THE INVENTION 

The invention discloses a method for predicting a fracture by analyzing at least one 
bone structure parameter. The method comprises: obtaining an image of a part of 
skeleton of a patient; locating at least one region of interest on the image of the 
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patient; extracting image data from the image of the patient; deriving at least one bone 
structure parameter from the image data of the patient; and predicting a fracture with 
the bone structure parameter of the patient. The bone structure parameter includes, 
but not limited to, bone micro-structure parameters and bone macro-structure 
parameters. 

[001 3] In certain aspects, described herein are methods of diagnosing, monitoring 
and/or predicting bone or articular disease (e.g., the risk of fracture) in a subject, the 
method comprises the steps of: determining one or more micro-structural parameters, 
and/or one or more macro-structure parameters, possibly with other bone parameters, 
of a bone or a joint in the subject; and combining at least two of the parameters to 
predict the risk of bone or articular disease. The micro-structural and macro-structure 
parameters may be, for example, one or more of the measurements/parameters shown 
in Tables 1 and 2. In certain embodiments, one or more micro-structural parameters 
and one or more macro-structural parameters are combined. In other embodiments, 
one or more micro-structural parameters and one or more other bone parameters are 
combined. In further embodiments, one or more macro-structure parameters and one 
or more other parameters are combined. In still further embodiments, one or more 
macro-structural parameters, one or more micro-structural parameters and one or 
more other bone parameters are combined. 

[0014] In any of the methods described herein, the comparing may comprise 
univariate, bivariate and/or multivariate statistical analysis of one or more of the 
parameters, including at least one bone structure parameter. In certain embodiments, 
the methods may further comprise comparing the parameters to data derived from a 
reference database of known disease parameters. 

[0015] In any of the methods described herein, the parameters are determined 
from an image obtained from the subject. In certain embodiments, the image 
comprises one or more regions of bone (e.g., patella, femur, tibia, fibula, pelvis, spine, 
etc). The image may be automatically or manually divided into two or more regions 
of interest. Furthermore, in any of the methods described herein, the image may be, 
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for example, an x-ray image, a CT scan, an MRI or the like and optionally includes 
one or more calibration phantoms. 



[001 6] In any of the methods described herein, the predicting includes performing 
univariate, bivariate or multivariate statistical analysis of the analyzed data and 
referencing the statistical analysis values to a fracture risk model. Fracture risk 
models can comprise, for example, data derived from a reference database of known 
fracture loads with their corresponding values of macro-anatomical, micro -anatomical 
parameters, and/or clinical risk factors. 

[001 7] In another aspect, the invention includes a method of determining the effect 
of a candidate agent on a subject's prognosis for musculoskeletal disease comprising: 
predicting a first risk of musculoskeletal disease in the subject according to any of the 
predictive methods described herein; administering a candidate agent to the subject; 
predicting a second risk of the musculoskeletal disease in the subject according to any 
of the predictive methods described herein; and comparing the first and second risks, 
thereby determining the effect of the candidate on the subject's prognosis for the 
disease. In any of these methods, the candidate agent can be administered to the 
subject in any modality, for example, by injection (intramuscular, subcutaneous, 
intravenous), by oral administration (e.g., ingestion), topical administration, mucosal 
administration or the like. Furthermore, the candidate agent may be a small molecule, 
a pharmaceutical, a biopharmaceutical, an agropharmaceuticals and/or combinations 
thereof. It is important to note that an effect on a subject's prognosis for 
musculoskeletal disease can occur in agents intended to have an effect, such as a 
therapeutic effect, on musculoskeletal disease as well as agents intended to primarily 
effect other tissues in the body but which have a secondary, or tangential, effect on 
musculoskeletal disease.. Further, the agent can be evaluated for the ability to effect 
diseases such as the risk of bone fracture (e.g., osteoporotic fracture). 

[0018] In other aspects, the invention includes a kit that is provided for aiding in 
the prediction of musculoskeletal disease (e.g., fracture risk). The kit typically 
comprises a software program that uses information obtained from an image to predict 
the risk or disease (e.g., fracture). The kit can also include a database of 
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measurements for comparison purposes. Additionally, the kit can include a subset of 
a database of measurements for comparisons. 

[001 9] In any of these methods, systems or kits, additional steps can be provided. 
Such additional steps include, for example, enhancing image data. 

[0020] Suitable subjects for these steps include for example mammals, humans 
and horses. Suitable anatomical regions of subjects include, for example, dental, 
spine, hip, knee and bone core x-rays. 

[0021] A variety of systems can be employed to practice the inventions. Typically 
at least one of the steps of any of the methods is performed on a first computer. 
Although, it is possible to have an arrangement where at least one of the steps of the 
method is performed on a first computer and at least one of the steps of the method is 
performed on a second computer. In this scenario the first computer and the second 
computer are typically connected. Suitable connections include, for example, a peer 
to peer network, direct link, intranet, and internet. 

[0022] It is important to note that any or all of the steps of the inventions 
disclosed can be repeated one or more times in series or in parallel with or without the 
repetition of other steps in the various methods. This includes, for example repeating 
the step of locating a region of interest, or obtaining image data. 

[0023] Data can also be converted from 2D to 3D to 4D and back; or from 2D to 
4D. Data conversion can occur at multiple points of processing the information. For 
example, data conversion can occur before or after pattern evaluation and/or analysis. 

[0024] Any data obtained, extracted or generated under any of the methods can be 
compared to a database, a subset of a database, or data previously obtained, extracted 
or generated from the subject. For example, known fracture load can be determined 
for a variety of subjects and some or all of this database can be used to predict 
fracture risk by correlating one or more micro-structural parameters or macro- 
structural parameters (Tables 1 and 2) with data from a reference database of fracture 
load for age, sex, race, height and weight matched individuals. 
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[0025] In any of the methods described herein, the analysis can comprise using 
one or more computer programs (or units). Additionally, the analysis can comprise 
identifying one or more regions of interest (ROI) in the image, either prior to, 
concurrently or after analyzing the image, e.g. for information on bone structure. 
Bone structural information can be, for example, one or more of the parameters 
shown in Table 1 and Table 2. The various analyses can be performed concurrently 
or in series. Further, when using two or more indices each of the indices can be 
weighted equally or differently, or combinations thereof where more than two indices 
are employed. Additionally, any of these methods can also include analyzing the 
image for bone structure information using any of the methods described herein. 

[0026] These and other embodiments of the subject invention will readily occur to 
those of skill in the art in light of the disclosure herein. 

BRIEF DESCRIPTION OF THE DRAWINGS 

[0027] Fig. 1 shows a flowchart of a method for collecting quantitative and/or 
qualitative data according to one embodiment of the present invention.. 

[0028] Fig.2 depicts exemplary regions of interest (ROIs), as analyzed in Example 
1 of the present invention. 

[0029] Figs.3A and 3B are graphs depicting correlation of 2D and 3D 
measurements according to one embodiment of the present invention. Fig. 3 A depicts 
correlation of 2D and 3D trabecular spacing. Fig. 3B depicts-correlation of 2D 
trabecular perimeter/trabecular area with 3D bone surface/bone volume. 

[0030] Figs.4A and 4B are graphs depicting correlation of 2D measurments with 
fracture load measurements according to one embodiment of the present invention. 
Fig. 4A depicts 2D trabecular perimeter/trabecular area v. fracture load. Fig. 4B 
depicts 2D trabecular separation vs. fracture load. 

[0031] Figs. 5 A and 5B are graphs depicting correlation of femoral neck DXA 
bone mineral density (BMD) and fracture load (Fig. 5A) and correlation of predicted 
fracture load and actual fracture load (Fig. 5B) according to a conventional method. 
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[0032] Figs.6A to 6D depict sliding window analysis maps for two different 
femora. The top maps (Figs. 6A and 6B) depict area ratio analysis. The bottom maps 
(Figs. 6C and 6D) depict trabecular perimeter analysis. Black lines show fracture 
lines from post-fracture x-rays for proximal and distal fragments. Red lines show 
results of watershed analysis of parameter maps. Color scale ranges from blue (low 
values) to red (high values). 

[0033] Figs.7A and 7B illustrate exemplary steps for predicting fracture risk via 
individualized fracture risk index (IFRI) according to one embodiment of the present 
invention. 

[0034] Fig. 8 shows an exemplary bone structure parameter map of a proximal 
femur. Color scale ranges from blue (low values) to red (high values). Regions are 
separated along low values ("valleys") using watershed segmentation. 

[0035] Fig.9 is an image of a femur and shows an approximation of the femoral 
axes by two linear segments in the neck (solid blue line) and shaft (solid green line). 
Also shown is a hyberbolic curve fitted to the intertrochanteric region. White arrows 
show loading simulating side impact fall. Three examples of cross-sectional lines are 
also shown. 

[0036] Fig. 10 depicts trochanteric and femoral neck fracture paths (red) that may 
be constructed by calculating the distance of segments (yellow) to cross-sectional 
lines (black).. 

[0037] Fig. 1 1 depicts definition of a region of interest (ROI) along the predicted 
fracture path using a region growing technique. This region of interest is used for a 
structural analysis of the trabecular bone. Contact points between the trabecular ROI 
and the cortical bone determine the area for cortical bone measurements. 

[0038] Fig. 12 is a block diagram of a computer program for predicting fracture 
risk according to one embodiment of the present invention. 
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DETAILED DESCRIPTION OF EMBODIMENTS 

[0039] The following description is presented to enable any person skilled in the 
art to make and use the invention. Various modifications to the embodiments 
described will be readily apparent to those skilled in the art, and the generic principles 
defined herein can be applied to other embodiments and applications without 
departing from the spirit and scope of the present invention as defined by the 
appended claims. Thus, the present invention is not intended to be limited to the 
embodiments shown, but is to be accorded the widest scope consistent with the 
principles and features disclosed herein. To the extent necessary to achieve a 
complete understanding of the invention disclosed, the specification and drawings of 
all issued patents, patent publications, and patent applications cited in this application 
are incorporated herein by reference. 

[0040] The practice of the present invention employs, unless otherwise indicated, 
methods of imaging and image processing within the skill of the art. Currently 
available imaging methods are explained fully in the literature. See, e.g., WO 
02/22014; X-Ray Structure Determination: A Practical Guide, 2 nd Edition, editors 
Stout and Jensen, 1989, John Wiley & Sons, publisher; Body CT: A Practical 
Approach, editor Slone, 1999, McGraw-Hill publisher; 77k? Essential Physics of 
Medical Imaging, editors Bushberg, Seibert, Leidholdt Jr & Boone, 2002, Lippincott, 
Williams & Wilkins; X-ray Diagnosis: A Physician } s Approach, editor Lam, 1998 
Springer- Verlag, publisher; Dental Radiology: Understanding the X-Ray Image, 
editor Laetitia Brocklebank 1997, Oxford University Press publisher; Digital Image 
Processing, editor Kenneth R. Castleman, 1996, Prentice Hall, publisher; The Image 
Processing Handbook, editor John C. Russ, 3 rd Edition, 1998, CRC Press; and Active 
Contours: The Application of Techniques from Graphics, Vision, Control Theory and 
Statistics to Visual Tracking of Shapes in Motion, Editors Andrew Blake, Michael 
Isard, 1999 Springer Verlag. As will be appreciated by those of skill in the art, as the 
field of imaging continues to advance, methods of imaging currently employed can 
evolve over time. Thus, any imaging method or technique that is currently employed 
is appropriate for application of the teachings of this invention as well as techniques 
that can be developed in the future. A further detailed description of imaging 
methods is not provided in order to avoid obscuring the invention. 
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[0041] Fig. 1 shows a method for collecting quantitative and/or qualitative data 
according to one embodiment of the present application. Step 101 is used to locate a 
part of the body of a subject, for example in a human body, for study. The part of the 
body located for study is the region of anatomical interest (RAI). In locating a part of 
the body for study, a determination is made to, for example, take an image or a series 
of images of the body at a particular location, e.g. hip, dental, spine, etc. Images 
include, for example, conventional x-ray images, x-ray tomosynthesis, ultrasound 
(including A-scan, B-scan and C-scan), computed tomography (CT scan), magnetic 
resonance imaging (MRI), optical coherence tomography, single photon emission 
tomography (SPECT), and positron emission tomography, or such other imaging tools 
that a person of skill in the art would find useful in practicing the invention. Once the 
image is taken, one or more regions of interest (ROI) can be manually and/or 
automatically located within the image at step 103. A skilled artisan would appreciate 
that algorithms can be used to automatically place regions of interest in a particular 
image. For instance, Example 1 below describes automatic placement of ROIs in 
femurs. Image data is extracted from the image at step 105. Finally, quantitative 
and/or qualitative data is extracted from the image data at step 107. The quantitative 
and/or qualitative data extracted from the image include at least one measurement 
about bone structure, such as those shown in Tables 1 and 2. 

[0042] Each step of locating a part of the body for study 101, optionally locating a 
region of interest 103, obtaining image data 105, and deriving quantitative and/or 
qualitative data 107, can be repeated one or more times at step 102,104, 106, or 108, 
respectively, as desired. Image data can be optionally enhanced by applying image 
processing techniques, such as noise filtering or diffusion filtering, to facilitate further 
analysis. 



11 



WO 2006/034018 



PCT/US2005/033131 



Table 1 

Representative Parameters Measured with 
Quantitative and Qualitative Image Analysis Methods for Micro-structure 



Parameter 


Measurements 


Measurements on 
extracted micro- 
structures 


•Trabecular contrast 

•Standard deviation of background subtracted ROI 
•Coefficient of Variation of ROI (Standard deviation / mean) 
•(Trabecular equivalent thickness / Marrow equivalent thickness) 
•Hough transform 
•Trabecular area 

(Pixel count of extracted trabeculae) 
•Trabecular area / Total area 
•Trabecular perimeter 

(Count of trabecular pixels with marrow pixels in their 

neighborhood, proximity or vicinity) 
•Trabecular distance transform 

(For each trabecular pixel, calculation of distance to closest 

marrow pixel) 
•Marrow distance transform 

(For each marrow pixel, calculation of distance to closest 

trabecular pixel) 
•Trabecular distance transform regional maximal values (mean, 

min., max, std. Dev). 

(Describes thickness and thickness variation of trabeculae) 
•Marrow distance transform regional maximal values (mean, min., 

max, std. Dev) 
•Star volume 

(Mean volume of all the parts of an object which can be seen 
unobscured from a random point inside the object in all possible 
directions) 
•Trabecular Bone Pattern Factor 

(TBPf = (PI - P2) / (Al - A2 ) where PI and Al are the 
perimeter length and trabecular bone area before dilation and P2 
and A2 corresponding values after a single pixel dilation, 
measure of connectivity) 


Measurements on 
skeleton of 
extracted micro- 
structures 


•Connected skeleton count or Trees (T) 
•Node count (N) 
•Segment count (S) 
•Node-to-node segment count (NN) 
•Node-to-free-end segment count (NF) 
•Node-to-node segment length (NNL) 
•Node-to-free-end segment length (NFL) 
•Free-end-to-free-end segment length (FFL) 
•Node-to-node total struts length (NNTSL) 
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Parameter 


Measurements 




•Free-end-to-free-ends total struts length( FF.TSL) 

•Total struts length (TSL) 

•FF.TSL/ TSL 

•NN.TSL/TSL 

•Loop count (Lo) 

•Loop area 

•Mean distance transform values for each connected skeleton 
•Mean distance transform values for each segment (Tb.Th ) 
•Mean distance transform values for each node-to-node segment 
(Tb.Th.NN) 

•Mean distance transform values for each node-to-free-end segment 

(Tb.Th.NF) 
•Orientation (angle) of each segment 
•Angle between segments 

•Length-thickness ratios (NNL/TbTh.NN ) and (NFL/ Tb.Th.NF) 
•Intercormectivity index (ICI) ICI = (N * NN)/ ( T * (NF + 1) ) 


Measurements on 
gray level images of 
micro-structures 


•Standard deviation of background subtracted ROI 
•Coefficient of Variation of ROI (Standard deviation / mean) 
•Fractal dimension 
•Fourier spectral analysis 

(Mean transform coefficient absolute value and mean spatial 

first moment) 
•Predominant orientation of spatial energy spectrum 

Watershed segmentation is applied to gray level images. Statistics 
of watershed segments are: 
•Total area of segments 

•Number of segments normalized by total area of segments 
•Average area of segments 
•Standard deviation of segment area 
•Smallest segment area 
•Largest segment area 



[0043] All micro-structural measurements can be applied in a direction-sensitive 
fashion or only on selected structures. For example, they can be applied to selected 
structures that are oriented parallel or perpendicular to stress lines. The techniques 
can also be used to measure only horizontal or vertical structures. 

[0044] As will be appreciated by those of skill in the art, the parameters and 
measurements shown in the tables are provided for illustration purposes only. It will 
be apparent that the terms micro-structural parameters, micro-architecture, micro- 
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anatomic structure, micro-structural and trabecular architecture may be used 
interchangeably. Furthermore, the terms macro-structural parameters, macro- 
structure, macro-anatomic parameters, macro-anatomic structure, macro-anatomy, 
macro-architecture and bone geometry may be used interchangeably. In addition, 
other parameters and measurements, ratios, derived values or indices can be used to 
extract quantitative and/or qualitative information about the ROI without departing 
from the scope of the invention. Additionally, where multiple ROI or multiple 
derivatives of data are used, the parameter measured can be the same parameter or a 
different parameter without departing from the scope of the invention. Additionally, 
data from different ROIs can be combined or compared as desired. Additional 
measurements can be performed that are selected based on the anatomical structure to 
be studied as described below. For instance, biomechanical aspects of the joint can 
also be evaluated. For example, the product of the average trabecular-computed 
tomography number and the total cross-sectional area of the sub-capital, basicervical 
or intertrochanteric regions can be determined, as it has been shown to correlate 
highly with failure loads. See, e.g., Lotz et al. (1990) J. Bone Joint Surg. Am. 72:689- 
700; Courtney et al. (1995) J. Bone Joint Surg. Am. 77(3):387-395; Pinilla et al. 
(1996) Calcif Tissue Int. 58:231-235. 

[0045] Once the quantitative and/or qualitative data is extracted from the image, it 
can be manipulated to assess the severity of the disease and to determine disease 
staging (e.g., mild, moderate, severe or a numerical value or index). The information 
can also be used to monitor progression of the disease and/or the efficacy of any 
interventional steps that have been taken. Finally, the information can be used to 
predict the progression of the disease or to randomize patient groups in clinical trials. 

[0046] After an image of an RAI is taken, one or more regions of interest can be 
identified within the image at step 103. The ROI can take up the entire image, or 
nearly the entire image. Alternatively, more than one ROI can be identified in an 
image, as shown in Fig. 2. One or more of the ROI may overlap or abut. As will be 
appreciated by a person of skill in the art, the number of ROI identified in an image is 
not limited to the seven depicted in Fig. 2. As also will be appreciated by those of 
skill in the art, where multiple ROI are used, any or all of the ROI can be organized 
such that it does not overlap, it abuts without overlapping, it overlaps partially, it 
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overlaps completely (for example where a first ROI is located completely within a 
second identified ROI), and combinations thereof. Further the number of ROI per 
image can range from one (ROIi) to n (ROI n ) where n is the number of ROI to be 
analyzed. 

[0047] Bone structure analyses, possibly together with bone density, and/or 
biomechanical (e.g. derived using finite element modeling) analyses, can be applied 
within a region of predefined size and shape and position. This region of interest can 
also be referred to as a "window". Processing can be applied repeatedly by moving 
the window to different positions of the image. For example, a field of sampling 
points can be generated and the analysis performed at these points. The results of the 
analyses for each parameter can be stored in a matrix space, e.g., where its position 
corresponds to the position of the sampling point where the analysis occurred, thereby 
forming a map of the spatial distribution of the parameter (a parameter map). The 
sampling field can have regular intervals or irregular intervals with varying density 
across the image. The window can have variable size and shape, for example to 
account for different patient size or anatomy. 

[0048] In another embodiment, rather than a fixed ROI (e.g., Fig. 3), the image 
may be overlaid with a regular grid, for example, a region of interest of a fixed size 
(e.g., of any shape) may be placed at each grid node, and parameters are evaluated 
within the boundaries of the ROI at each position. This results in a value for each 
bone parameter at each grid node, which can be displayed in a color-coded map of the 
proximal femur for each parameter. 

[0049] The amount of overlap between the windows can be determined, for 
example, using the interval or density of the sampling points (and resolution of the 
parameter maps). Thus, the density of sampling points is set higher in regions where 
higher resolution is desired and set lower where moderate resolution is sufficient, in 
order to improve processing efficiency. The size and shape of the window would 
determine the local specificity of the parameter. Window size is preferably set such 
that it encloses most of the structure being measured. Oversized windows are 
generally avoided to help ensure that local specificity is not lost. 
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[0050] The shape of the window can be varied to have the same orientation and/or 
geometry of the local structure being measured to minimize the amount of structure 
clipping and to maximize local specificity. Thus, both 2D and/or 3D windows can be 
used, as well as combinations thereof, depending on the nature of the image and data 
to be acquired. 

[0051] In another embodiment, bone structure analysis, possibly together with 
bone density and/or biomechanical (e.g. derived using finite element modeling) 
analyses, can be applied within a region of predefined size and shape and position. 
The region is generally selected to include most, or all, of the anatomic region under 
investigation and, preferably, the parameters can be assessed on a pixel-by-pixel basis 
(e.g., in the case of 2D or 3D images) or a voxel-by- voxel basis in the case of cross- 
sectional or volumetric images (e.g., 3D images obtained using MR and/or CT). 
Alternatively, the analysis can be applied to clusters of pixels or voxels wherein the 
size of the clusters is typically selected to represent a compromise between spatial 
resolution and processing speed. Each type of analysis can yield a parameter map. 

[0052] Parameter maps can be based on measurement of one or more parameters 
in the image or window; however, parameter maps can also be derived using 
statistical methods. In one embodiment, such statistical comparisons can include 
comparison of data to a reference population, e.g. using a z-score or a T-score. Thus, 
parameter maps can include a display of z-scores or T- scores. 

[0053] Additional measurements relating to the site to be measured can also be 
taken. For example, measurements can be directed to dental, spine, hip, knee or bone 
cores. Examples of suitable site specific measurements are shown in Table 2. 
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Table 2 

Common and site specific measurements of bone macro-structure parameters 



Measurements on 
macro-structures 
common to dental, 
spine, hip, knee or bone 
cores images 


The following parameters are derived from the extracted 
macro-structures : 

•Calibrated density of extracted structures 
•Calibrated density of background 
•Average intensity of extracted structures 
•Average intensity of background (area other than extracted 
structures) 

•Structural contrast (average intensity of extracted structures / 

average intensity of background ) 
•Calibrated structural contrast (calibrated density extracted 

structures / calibrated density of background) 
•Total area of extracted structures 

•Bone pattern factor; measures concavity and convexity of 
structures 

•Average length of structures (units of connected segments) 
•Maximum length of structures 
•Average thickness of structures 
•Maximum thickness of structures 
•Regional maximum thickness of structures 
•Standard deviation of thickness along structures 
•Average orientation angle of structure segmentss 
•Structure segment tortuosity; a measure of straightness 
•Structure segment solidity; another measure of straightness 


Parameters specific to 
hip images 


•Shaft angle 
•Neck angle 

•Average and minimum diameter of femur neck 

w J. 1 Y vl v4AAVA 11111 A J, A J. X Vill 1 V4-A WJLJLJL V IVl J.V111VU. IIVVIV 

•Hip axis length 

•CCD (caput-collum-diaphysis) angle 

•Width of trochanteric region 

•Largest cross-section of femur head 

•Standard deviation of cortical bone thickness within ROI 

•Minimum, maximum, mean and median thickness of cortical 

bone within ROI 
•Hip j oint space width 


Parameters specific to 
snine images 


•Superior endplate cortical thickness (anterior, center, 

•Inferior endplate cortical thickness (anterior, center, posterior) 
•Anterior vertebral wall cortical thickness (superior, center, 
inferior) 

•Posterior vertebral wall cortical thickness (superior, center, 
inferior) 

•Superior aspect of pedicle cortical thickness 
•inferior aspect of pedicle cortical thickness 
•Vertebral height (anterior, center, posterior) 
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•Vertebral diameter (superior, center, inferior), 




•Pedicle thickness (supero-inferior direction). 




•Maximum vertebral height 




•Minimum vertebral height 




•Average vertebral height 




•Anterior vertebral height 




•Medial vertebral height 




•Posterior vertebral height 




•Maximum inter-vertebral height 




•Minimum inter-vertebral height 




•Average inter-vertebral height 


Parameters specific to 


•Average medial joint space width 


knee images 


•Minimum medial joint space width 




•Maximum medial joint space width 




•Average lateral joint space width 




•Minimum lateral joint space width 




•Maximum lateral joint space width 



[0054] As will be appreciated by those of skill in the art, measurement and image 
processing techniques are adaptable to be applicable to both micro-architecture and 
macro-anatomical structures. Examples of these measurements are shown in Table 2, 

[0055] As noted above, analysis can also include one or more additional 
techniques,, for example, Hough transform, mean pixel intensity analysis, variance of 
pixel intensity analysis, soft tissue analysis and the like. See, e.g., co-owned 
International Application WO 02/30283. 

[0056] Calibrated density typically refers to the measurement of intensity values of 
features in images converted to its actual material density or expressed as the density 
of a reference material whose density is known. The reference material can be metal, 
polymer, plastics, bone, cartilage, etc., and can be part of the object being imaged or a 
calibration phantom placed in the imaging field of view during image acquisition. 

[0057] Extracted structures typically refer to simplified or amplified 
representations of features derived from images. Bone structure parameters include, 
for example, micro-structure parameters and macro-structure parameters. Micro- 
structure parameters could be, for example , the measurements in Table 1 . Macro- 
structure parameters could be, for example, the parameters in Table 2. An example 
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would be binary images of trabecular patterns generated by background subtraction 
and thresholding. Another example would be binary images of cortical bone 
generated by applying an edge filter and thresholding. The binary images can be 
superimposed on gray level images to generate gray level patterns of structure of 
interest. 

[0058] Distance transform typically refers to an operation applied on binary 
images where maps representing distances of each 0 pixel to the nearest 1 pixel are 
generated. Distances can be calculated by the Euclidian magnitude, city-block 
distance, La Place distance or chessboard distance. 

[0059] Distance transform of extracted structures typically refers to distance 
transform operation applied to the binary images of extracted structures, such as those 
discussed above with respect to calibrated density. 

[0060] Skeleton of extracted structures typically refers to a binary image of 1 pixel 
wide patterns, representing the centerline of extracted structures. It is generated by 
applying a skeletonization or medial transform operation, by mathematical 
morphology or other methods, on an image of extracted structures. 

[0061] Skeleton segments typically are derived from skeleton of extracted 
structures by performing pixel neighborhood analysis on each skeleton pixel. This 
analysis classifies each skeleton pixel as a node pixel or a skeleton segment pixel. A 
node pixel has more than 2 pixels in its 8 -neighborhood. A skeleton segment is a 
chain of skeleton segment pixels continuously 8 -connected. Two skeleton segments 
are separated by at least one node pixel. 

[0062] Watershed segmentation as it is commonly known to a person of skill in the 
art, typically is applied to gray level images to characterize gray level continuity of a 
structure of interest. The statistics of dimensions of segments generated by the 
process are, for example, those listed in Table 1 above. As will be appreciated by 
those of skill in the art, however, other processes can be used without departing from 
the scope of the invention. As described in the Examples, watershed transformation 
may be applied as follows. The image (or its negative, depending on whether peaks 
or valleys are to be located) is considered as a topographic relief, in which higher 
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intensities correspond to higher topographic heights. This relief can be divided 
(segmented) into catchment basins, one for each local minimum of the image, where a 
catchment basin is defined as the area in which a raindrop would flow to the 
corresponding minimum. The lines that separate catchment basins from each other 
are the watersheds. 

[0063] At step 109, the extracted image data obtained at step 107 can be converted 
to a 2D pattern, a 3D pattern or a 4D pattern, for example including velocity or time, 
to facilitate data analyses. Following conversion to 2D, 3D or 4D pattern the images 
are evaluated for patterns at step 111. Additionally images can be converted from 2D 
to 3D, or from 3D to 4D, if desired according to step 110. Persons of skill in the art 
will appreciate that similar conversions can occur between 2D and 4D in this process 
or any process illustrated in this invention. 

[0064] As will be appreciated by those of skill in the art, the conversion step is 
optional and the process can proceed directly from extracting image data from the 
ROI at step 107 to evaluating the data pattern at step 111. Evaluating the data for 
patterns, includes, for example, performing the measurements described in Table 1 or 
Table 2, above. 

[0065] Additionally, the steps of locating the region of interest, obtaining image 
data, and evaluating patterns can be performed once or a plurality of times, 
respectively at any stage of the process. For example, following an evaluation of 
patterns at step 111, additional image data can be obtained according to step 1 14, or 
another region of interest can be located according to step 112. These steps can be 
repeated as often as desired, in any combination desirable to achieve the data analysis 
desired. 

[0066] An alternative process includes the step of enhancing image data prior to 
converting an image or image data to a 2D, 3D, or 4D pattern. The process of 
enhancing image data, can be repeated if desired. In still further embodiments, the 
step of enhancing image data may occur after converting an image or image data to a 
2D, 3D, or 4D pattern. Again, the process of enhancing image data , can be repeated 
if desired. 
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[0067] Furthermore, in certain embodiments, after locating a part of the body for 
study and imaging, the image is then converted to a 2D pattern, 3D pattern or 4D 
pattern- The region of interest is optionally located within the image after optional 
conversion to a 2D, 3D and/or 4D image and data is then extracted. Patterns are then 
evaluated in the extracted image data 

[0068] Some or all the processes can be repeated one or more times as desired. 
For example, locating a part of the body for study, locating a region of interest, 
obtaining image data, and evaluating patterns, can be repeated one or more times if 
desired, respectively. For example, following an evaluation of patterns, additional 
image data can be obtained, or another region of interest can be located and/or 
another portion of the body can be located for study. These steps can be repeated as 
often as desired, in any combination desirable to achieve the data analysis desired. 

[0069] Image data may also be enhanced The step of enhancing image data may 
occur prior to conversion, prior to locating a region of interest, prior to obtaining 
image data, or prior to evaluating patterns. 

[0070] The method also comprises obtaining an image of a bone or a joint, 
optionally converting the image to a two-dimensional or three-dimensional or four- 
dimensional pattern, and evaluating the amount or the degree of normal, diseased or 
abnormal tissue or the degree of degeneration in a region or a volume of interest using 
one or more of the parameters specified in Table 1 and/or Table 2. By performing 
this method at an initial time Ti, information can be derived that is useful for 
diagnosing one or more conditions or for staging, or determining, the severity of a 
condition. This information can also be useful for determining the prognosis of a 
patient, for example with osteoporosis or arthritis. By performing this method at an 
initial time Ti, and a later time T2, the change, for example in a region or volume of 
interest, can be determined which then facilitates the evaluation of appropriate steps 
to take for treatment. Moreover, if the subject is already receiving therapy or if 
therapy is initiated after time Ti, it is possible to monitor the efficacy of treatment. 
By performing the method at subsequent times, T2-T n . additional data can be acquired 
that facilitate predicting the progression of the disease as well as the efficacy of any 
interventional steps that have been taken. As will be appreciated by those of skill in 
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the art, subsequent measurements can be taken at regular time intervals or irregular 
time intervals, or combinations thereof. For example, it caln be desirable to perform 
the analysis at Ti with an initial follow-up, T 2 , measurement taken one month later. 
The pattern of one month follow-up measurements could be performed for a year (12 
one-month intervals) with subsequent follow-ups performed at 6 month intervals and 
then 12 month intervals. Alternatively, as an example, three initial measurements 
could be at one month, followed by a single six month follow up which is then 
followed again by one or more one month follow-ups prior to commencing 12 month 
follow ups. The combinations of regular and irregular intervals are endless, and are 
not discussed further to avoid obscuring the invention. 

[0071] Moreover, one or more of the bone structure parameters listed in Tables 1 
and 2, and possibly one or more parameters, can be measured. The measurements can 
be analyzed separately or the data can be combined, for example using statistical 
methods such as linear regression modeling or correlation. Actual and predicted 
measurements can be compared and correlated. See, also, Examples described later. 

[0072] The method for predicting future fracture in a subject can be fully 
automated such that the measurements of one or more of the bone structure 
parameters specified in Tables 1 and 2, and possibly one or more other parameters, 
are done automatically without intervention. As will be appreciated by those of skill 
in the art, the fully automated analysis is, for example, possible with one or more of 
the steps involved in predicting future fracture, including, sliding window ROI 
analysis of such bone parameter(s) to generate bone parameter maps; watershed 
segmentation of parameter maps to identify possible or likely fracture lines; local 
structure analysis (e.g., placement of ROI along predicted fracture line and analysis of 
traceular and cortical bone parameters); combining multiple bone parameters; and/or 
calculations such as multivariate regressions. This process may also include, for 
example, seed growing, thresholding, atlas and model bas^d segmentation methods, 
live wire approaches, active and/or deformable contour approaches, contour tracking, 
texture based segmentation methods, rigid and non-rigid surface or volume 
registration, for example based on mutual information or other similarity measures. 
One skilled in the art will readily recognize other techniques and methods for folly 
automated assessment of the parameters and measurements described herein. 
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[0073] Alternatively, the method of predicting future fractures in a subject can be 
semi-automated such that the measurements of one or more of the parameters, 
including at least one bone structure parameter, are performed semi-automatically, 
i.e., with intervention. The semi-automatic assessment allows for human interaction 
and, for example, quality control, and utilizing the measurement of such parameter(s) 
to diagnose, stage, prognosticate or monitor a disease or to monitor a therapy. The 
semi-automated measurement is, for example, possible with image processing 
techniques such as segmentation and registration. This can include seed growing, 
thresholding, atlas and model based segmentation methods, live wire approaches, 
active and/or deformable contour approaches, contour tracking, texture based 
segmentation methods, rigid and non-rigid surface or volume registration, for example 
base on mutual information or other similarity measures. One skilled in the art will 
readily recognize other techniques and methods for semi-automated assessment of 
such parameters. 

[0074] Following the step of deriving quantitative and/or qualitative image data, 
one or more candidate agents can be administered to the patient . The candidate agent 
can be any agent the effects of which are to be studied. Agents can include any 
substance administered or ingested by a subject, for example, molecules, 
pharmaceuticals, biopharmaceuticals, agropharmaceuticals, or combinations thereof, 
including cocktails, that are thought to affect the quantitative and/or qualitative 
parameters that can be measured in a region of interest. These agents are not limited 
to those intended to treat disease that affects the musculoskeletal system but this 
invention is intended to embrace any and all agents regardless of the intended 
treatment site. Thus, appropriate agents are any agents whereby an effect can be 
detected via imaging. The steps of locating a region of interest, obtaining image data, 
obtaining such quantitative and/or qualitative data from image data, and administering 
a candidate agent, can be repeated one or more times as desired, respectively. Image 
data may be enhanced as often as desired. 

[0075] Furthermore, an image may be taken prior to administering the candidate 
agent. However, as will be appreciated by those of skill in the art, it is not always 
possible to have an image prior to administering the candidate agent. In those 
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situations, progress is determined over time by evaluating the change in parameters 
from extracted image to extracted image. 



[0076] The derived quantitative and/or qualitative information can be compared to 
an image taken at Tl, or any other time, if such image is available. Again, the steps 
of deriving information and/or enhancing data can be repeated, as desired. 

[0077] In addition, following the step of extracting image data from the ROI, the 
image can be transmitted Transmission can be to another computer in the network or 
via the World Wide Web to another network. Following the step of transmitting the 
image, the image is converted to a pattern of normal and diseased tissue. Normal 
tissue includes the undamaged tissue located in the body part selected for study. 
Diseased tissue includes damaged tissue located in the body part selected for study. 
Diseased tissue can also include, or refer to, a lack of normal tissue in the body part 
selected for study. For example, damaged or missing bone would be considered 
diseased tissue. Once the image is converted, it may be analyzed. 

[0078] The step of transmitting the image is optional. As will be appreciated by 
those of skill in the art, the image can also be analyzed prior to converting the image 
to a pattern of normal and diseased. 

[0079] As previously described, some or all the processes can be repeated one or 
more times as desired. For example, locating a region of interest, obtaining image 
data, enhancing image data, transmitting an image, converting the image to a pattern 
of normal and diseased, analyzing the converted image, can be repeated one or more 
times if desired, respectively. 

[0080] Two or more devices may be connected. Either the first or second device 
can develop a degeneration pattern from an image of a region of interest. Similarly, 
either device can house a database for generating additional patterns or measurements. 
The first and second devices can communicate with each other in the process of 
analyzing an image, developing a degeneration pattern from a region of interest in the 
image, creating a dataset of patterns or measurements or comparing the degeneration 
pattern to a database of patterns or measurements. However, all processes can be 
performed on one or more devices, as desired or necessary. 
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[0081] In this method the electronically generated, or digitized image or portions 
of the image can be electronically transferred from a transferring device to a receiving 
device located distant from the transferring device; receiving the transferred image at 
the distant location; converting the transferred image to a pattern of normal or 
diseased or abnormal tissue using one or more of the parameters specified in Table 1 
or Table 2; and optionally transmitting the pattern to a site for analysis. As will be 
appreciated by those of skill in the art, the transferring device and receiving device 
can be located within the same room or the same building. The devices can be on a 
peer-to-peer network, or an intranet. Alternatively, the devices can be separated by 
large distances and the information can be transferred by any suitable means of data 
transfer, including http and ftp protocols. 

[0082] Alternatively, the method can comprise electronically transferring an 
electronically-generated image or portions of an image of a bone or a joint from a 
transferring device to a receiving device located distant from the transferring device; 
receiving the transferred image at the distant location; converting the transferred 
image to a degeneration pattern or a pattern of normal or diseased or abnormal tissue 
using one or more of the parameters specified in Table 1 or Table 2; and optionally 
transmitting the degeneration pattern or the pattern of normal or diseased or abnormal 
tissue to a site for analysis. 

[0083] Thus, the invention described herein includes methods and systems for 
prognosis of fracture risk. (See, also, Examples). 

[0084] In order to make more accurate prognoses, it may be desirable in certain 
instances to compare data obtained from a subject to a reference database. For 
example, when predicting fracture risk, it may be useful to compile data of actual 
(known) fracture load in a variety of samples and store the results based on clinical 
risk factors such as age, sex and weight (or other characteristics) of the subject from 
which the sample is obtained. The images of these samples are analyzed to obtain 
parameters shown in Tables 1 and 2, and possibly one or more other parameters,. A 
fracture risk model correlated with fracture load may be developed using univariate, 
bivariate and/or multivariate statistical analysis of these parameters and is stored in 
this database. A fracture risk model may include information that is used to estimate 
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fracture risk from parameters shown in Tables 1 and 2, and possibly one or more other 
parameters,. An example of a fracture risk model is the coefficients of a multivariate 
linear model derived from multivariate linear regression of these parameters (Tables 
1, 2, age, sex, weight, etc.) with fracture load. A person skilled in the art will 
appreciate that fracture risk models can be derived using other methods such as 
artificial neural networks and be represented by other forms such as the coefficients of 
artificial neural networks. Patient fracture risk can then be determined from 
measurements obtain from bone images by referencing to this database. 

[0085] In conventional methods of determining actual fracture load, cross- 
sectional images may be taken throughout testing to determine at what load force a 
fracture might occur. 

[0086] The analysis techniques described herein can then be applied to a subject 
and the risk of fracture (or other disease) could be predicted using one or more of the 
parameters described herein. The prognostication methods described herein are more 
accurate than conventional methods for predicting fracture risk. Fig. 5 A is a graph 
depicting conventional linear regression analysis of DXA bone mineral density 
correlated to fracture load. Correlations of individual parameters to fracture load are 
comparable to DXA. However, when multiple structural parameters are combined, 
the prediction of load at which fracture will occur is more accurate. Thus, the 
analyses of images as described herein can be used to accurately predict 
musculoskeletal disease such as fracture risk. 

[0087] Another aspect of the present invention is a kit for aiding in predicting 
fracture risk in a subject, which kit comprises a software program, which when 
installed and executed on a computer creates a bone parameter map (e.g., using one or 
more of the parameters specified in Tables 1 and 2, and possibly one or more other 
parameters) presented in a standard graphics format and produces a computer readout. 
The kit can further include software for (1) identifying likely fracture lines (e.g., by 
watershed segmentation); (2) placing one or more ROI along predicted fracture 
line(s); (3) analyzing one or more bone parameters along predicted fracture lines; 
and/or (4) combining multiple bone parameters and calculating fracture load. 
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[0088] The kit can further include one or more databases of measurements for use 
in calibrating or diagnosing the subject. One or more databases can be provided to 
enable the user to compare the results achieved for a specific subject against, for 
example, a wide variety of subjects, or a small subset of subjects having 
characteristics similar to the subject being studied. 

[0089] A system is provided that includes (a) a device for electronically 
transferring an image, a parameter map, an analyzed parameter map, etc., to a 
receiving device located distant from the transferring device; (b) a device for 
receiving the image or map at the remote location; (c) a database accessible at the 
remote location for generating additional patterns or measurements for the bone or the 
joint of a subject wherein the database includes a collection of subject patterns or 
data, for example of human bones or joints, which patterns or data are organized and 
can be accessed by reference to characteristics such as type of joint, gender, age, 
height, weight, bone size, type of movement, and distance of movement; and (d) 
optionally a device for transmitting the correlated pattern back to the source of the 
degeneration pattern or pattern of normal, diseased or abnormal tissue. 

[0090] Thus, the methods and systems described herein may make use of 
collections of data sets of measurement values, for example measurements of bone 
structure, probably with other measurements from images (e.g., x-ray images). 
Records can be formulated in spreadsheet-like format, for example including data 
attributes such as date of image (x-ray), patient age, sex, weight, current medications, 
geographic location, etc. The database formulations can further comprise the 
calculation of derived or calculated data points from one or more acquired data points, 
typically using the parameters listed in Tables 1 and 2 or combinations thereof. A 
variety of derived data points can be useful in providing information about individuals 
or groups during subsequent database manipulation, and are therefore typically 
included during database formulation. Derived data points include, but are not limited 
to the following: (1) maximum value of a selected bone structure parameter, 
determined for a selected region of bone or joint or in multiple samples from the same 
or different subjects; (2) minimum value of a selected bone structure parameter, 
determined for a selected region of bone or joint or in multiple samples from the same 
or different subjects; (3) mean value of a selected bone structure parameter, 
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determined for a selected region of bone or joint or in multiple samples from the same 
or different subjects; (4) the number of measurements that are abnormally high or 
low, determined by comparing a given measurement data point with a selected value; 
and the like. Other derived data points include, but are not limited to the following: 
(1) maximum value of bone mineral density, determined for a selected region of bone 
or in multiple samples from the same or different subjects; (2) minimum value of 
bone mineral density, determined for a selected region of bone or in multiple samples 
from the same or different subjects; (3) mean value of bone mineral density, 
determined for a selected region of bone or in multiple samples from the same or 
different subjects; (4) the number of bone mineral density measurements that are 
abnormally high or low, determined by comparing a given measurement data point 
with a selected value; and the like. Other derived data points will be apparent to 
persons of ordinary skill in the art in light of the teachings of the present specification. 
The amount of available data and data derived from (or arrived at through analysis of) 
the original data provides an unprecedented amount of information that is very 
relevant to management of musculoskeletal-related diseases such as osteoporosis or 
arthritis. For example, by examining subjects over time, the efficacy of medications 
can be assessed. 

[0091] Measurements and derived data points are collected and calculated, 
respectively, and can be associated with one or more data attributes to form a 
database. 

[0092] Data attributes can be automatically input with the electronic image and 
can include, for example, chronological information (e.g., DATE and TIME). Other 
such attributes can include, but are not limited to, the type of imager used, scanning 
information, digitizing information and the like. Alternatively, data attributes can be 
input by the subject and/or operator, for example subject identifiers, i.e., 
characteristics associated with a particular subject. These identifiers include but are 
not limited to the following: (1) a subject code (e.g., a numeric or alpha-numeric 
sequence); (2) demographic information such as race, gender and age; (3) physical 
characteristics such as weight, height and body mass index (BMI); (4) selected 
aspects of the subject's medical history (e.g., disease states or conditions, etc.); and 
(5) disease-associated characteristics such as the type of bone disorder, if any; the 
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type of medication used by the subject. In the practice of the present invention, each 
data point would typically be identified with the particular subject, as well as the 
demographic, etc. characteristic of that subject. 

[0093] Other data attributes will be apparent to persons of ordinary skill in the art 
in light of the teachings of the present specification. (See, also, WO 02/30283). 

[0094] Thus, data about bone structure information, possibly with bone mineral 
density information and/or articular information, is obtained from normal control 
subjects using the methods described herein. These databases are typically referred to 
as "reference databases" and can be used to aid analysis of any given subject's image, 
for example, by comparing the information obtained from the subject to the reference 
database. Generally, the information obtained from the normal control subjects will 
be averaged or otherwise statistically manipulated to provide a range of "normal" 
measurements. Suitable statistical manipulations and/or evaluations will be apparent 
to those of skill in the art in view of the teachings herein. The comparison of the 
subject's information to the reference database can be used to determine if the 
subject's bone information falls outside the normal range found in the reference 
database or is statistically significantly different from a normal control. 

[0095] Data obtained from images, as described above, can be manipulated, for 
example, using a variety of statistical analyses to produce useful information. 
Databases can be created or generated from the data collected for an individual, or for 
a group of individuals, over a defined period of time (e.g., days, months or years),, 
from derived data, and from data attributes. 

[0100] For example, data can be aggregated, sorted, selected, sifted, clustered and 
segregated by means of the attributes associated with the data points. A number of 
data mining software exist which can be used to perform the desired manipulations. 

[0101] Relationships in various data can be directly queried and/or the data 
analyzed by statistical methods to evaluate the information obtained from 
manipulating the database. 
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[0102] For example, a distribution curve can be established for a selected data set, 
and the mean, median and mode calculated therefor. Further, data spread 
characteristics, e.g., variability, quartiles, and standard deviations can be calculated. 

[0103] The nature of the relationship between any variables of interest can be 
examined by calculating correlation coefficients. Useful methods for doing so 
include, but are not limited to: Pearson Product Moment Correlation and Spearman 
Rank Correlation. Analysis of variance permits testing of differences among sample 
groups to determine whether a selected variable has a discernible effect on the 
parameter being measured. 

[01 04] Non-parametric tests can be used as a means of testing whether variations 
between empirical data and experimental expectancies are attributable to chance or to 
the variable or variables being examined. These include the Chi Square test, the Chi 
Square Goodness of Fit, the 2x2 Contingency Table, the Sign Test and the Phi 
Correlation Coefficient. Other tests include z-scores, T-scores or lifetime risk for 
arthritis, cartilage loss or osteoporotic fracture. 

[0105] There are numerous tools and analyses available in standard data mining 
software that can be applied to the analyses of the databases that can be created 
according to this invention. Such tools and analysis include, but are not limited to, 
cluster analysis, factor analysis, decision trees, neural networks, rule induction, data 
driven modeling, and data visualization. Some of the more complex methods of data 
mining techniques are used to discover relationships that are more empirical and data- 
driven, as opposed to theory driven, relationships. 

[01 06] Statistical significance can be readily determined by those of skill in the 
art. The use of reference databases in the analysis of images facilitates that diagnosis, 
treatment and monitoring of bone conditions such as osteoporosis. 

[0107] For a general discussion of statistical methods applied to data analysis, see 
Applied Statistics for Science and Industry, by A. Romano, 1977, Allyn and Bacon, 
publisher. 
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[01 08] The data is preferably stored and manipulated using one or more computer 
programs or computer systems. These systems will typically have data storage 
capability (e.g., disk drives, tape storage, optical disks, etc.)- Further, the computer 
systems can be networked or can be stand-alone systems. If networked, the computer 
system would be able to transfer data to any device connected to the networked 
computer system for example a medical doctor or medical care facility using standard 
e-mail software, a central database using database query and update software (e.g., a 
data warehouse of data points, derived data, and data attributes obtained from a large 
number of subjects). Alternatively, a user could access from a doctor's office or 
medical facility, using any computer system with Internet access, to review historical 
data that can be useful for determining treatment. 

[01 09] If the networked computer system includes a World Wide Web 
application, the application includes the executable code required to generate database 
language statements, for example, SQL statements. Such executables typically 
include embedded SQL statements. The application further includes a configuration 
file that contains pointers and addresses to the various software entities that are 
located on the database server in addition to the different external and internal 
databases that are accessed in response to a user request. The configuration file also 
directs requests for database server resources to the appropriate hardware, as can be 
necessary if the database server is distributed over two or more different computers. 

[0110] As a person of skill in the art will appreciate, one or more of the 
parameters specified in Table 1 and Table 2 can be used at an initial time point Ti to 
assess the severity of a bone disease such as osteoporosis. The patient can then serve 
as their own control at a later time point T2, when a subsequent measurement using 
one or more of the same parameters used at Ti is repeated. 

[0111] A variety of data comparisons can be made that will facilitate drug 
discovery, efficacy, dosing, and comparisons. For example, one or more of the 
parameters specified in Table 1 and Table 2 may be used to identify lead compounds 
during drug discovery. For example, different compounds can be tested in animal 
studies and the lead compounds with regard to the highest therapeutic efficacy and 
lowest toxicity, e.g. to the bone or the cartilage, can be identified. Similar studies can 
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be performed in human subjects, e.g., FDA phase I, II or III trials. Alternatively, or in 
addition, one or more of the parameters specified in Table 1 and Table; 2 can be used 
to establish optimal dosing of a new compound. It will be appreciated also that one or 
more of the parameters specified in Table 1 and Table 2 can be used to compare a 
new drug against one or more established drugs or a placebo. The patient can then 
serve as their own control at a later time point T2. 

EXAMPLES 

Example 1: Correlation of Micro-Structural and Macro-Structural 
Parameters to Fracture Load 

[0112] Using 1 5 fresh cadaveric femurs, the following analyses were performed 
to determine the correlation of various micro-structural and macro-structural 
parameters to fracture load, as determined by biomechanical testing. Parameters 
measured included one or more of the following 



Parameter Name 


Description 


Measurements on gray values of extracted structures 


Std. dev. of 
normalized ROI 


Normalized ROI is subtracted from the background using a difference of 
gaussian filter. The standard deviation reflects the "roughness" of the 
trabecular structures. 


Measurements on binarization of extracted structures 


Trab. Perimeter 


Total length of outline (perimeter) of extracted trabecular structures in a 
ROI. 


Trab. 

Perimeter/Trab. 
Area 


Trabecular perimeter normalized by area of extracted trabecular structures. 


Trab. 

Perimeter/Total 
Area 


Trabecular perimeter normalized by ROI area. 


Trabecular Bone 
Pattern Factor 


Change of perimeter per change of area. Measures concavity and convexity 
of structures. 


Trabecular Star 
Volume 


Estimated volume of trabecular structures by measuring distance of random 
points to boundaries of extracted structures. 


Marrow Space 
Length 


Mean length of skeletonized marrow space (background) region. 
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Parameter Name 


Description 


Mode Trab. 
Separation 


The mode of distance transform values of the marrow space (background) 
region. 


Std. Dev of 

Trabecular 

Separation 


The standard deviation of distance transform values of the marrow space 
(background) region. 


Trabecular 
Separation 


The mean of distance transform values of the marrow space region. 


Trabecular 
Thickness 


The mean of distance transform values along the skeleton (centerline) of 
extracted structures. 


Max.Trab. 
Thickness 


The maximum distance transform value of extracted structures in an ROI. 


Measurements on skeleton of extracted structures 


Trabecular 

Segment 

Thickness 


The mean of distance transform values along the segmented (by nodes) 
skeleton of extracted structures. 


Free-end Segment 
Thickness 


The mean of distance transform values along the free-end segments of the 
skeletonized structures. 


Node-Node 

Segment 

Thickness 


The mean of distance transform values along the node-node (inner) 
segments of the skeletonized structures. 


Number of Nodes 


Number of nodes (branching points) of skeletonized structures normalized 
by ROI area. 


Segment Number 


Number of skeleton segments normalized by ROI area. 


Free-end Segment 
Number 


Number of free-end skeleton segments normalized by ROI area. 


Segment 
Tortuosity 


Ratio of length of segments to distance between segment ends. 


Segment Solidity 


Ratio of length of segment to area of convex hull of the segment. 


Watershed segmentation is applied to normalized gray level images. Statistics of 
watershed segments are: 


Watershed 
Segment Area 


Average area of watershed segments. Measures the trabecular separation by 
area between structures. 


Watershed 
Segment Number 


Number of watershed segments normalized by ROI area. 


Std. dev. of 
Watershed Area 


Standard deviation of areas of watershed segments. Measures the 
homogeneity of trabecular separation 


Macro-anatomical and geometric parameters 
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Parameter Name 


Description 


Median Cortical 
Thickness 


Median of distance transform values measured along the centerline of 
extracted cortical bone structure. 


Maximum 
Cortical 
i juckucss 


Maximum of distance transform values measured along the centerline of 
extracted cortical bone structure. 


nip /\X1S .uengul 


Length of the femoral neck axis, extending from the bone edge at the base 
of trochanter to the bone edge at the inner pelvic brim (femoral head for 
cadaveric femur). 


Neck-shaft angle 


Angle between femoral neck axis and shaft axis. 


Head diameter 




Mean Neck Width 


Mean of distance transform on femoral neck axis between center of femoral 
head to intertrochanteric line. 


Minimum Neck 
Width 


Minimum distance transform value on femoral neck axis. 



[0113] Standardization of Hip radiographs : Density and magnification 
calibration on the x-ray radiographs was achieved using a calibration phantom. The 
reference orientation of the hip x-rays was the average orientation of the femoral 
shaft. 

[0114] Automatic Placement of Regions of Interest. Seven regions of interest 
were consistently and accurately placed based on the geometry and position of the 
proximal femur (Fig. 2). This was achieved by detecting femoral boundaries, 
estimating shaft and neck axes, and constructing the ROIs based on axes and 
boundary intercept points. This approach ensured that the size and shape of ROIs 
placed conformed to the scale and shape of the femur, and thus were consistent 
relative to anatomic features on the femur. 

[0115] Automatic Segmentation of the proximal femur: A global gray level 
thresholding using bi-modal histogram segmentation algorithm(s) was performed on 
the hip images and a binary image of the proximal femur was generated. Edge- 
detection analysis was also performed on the hip x-rays, including edge detection of 
the outline of the proximal femur that involved breaking edges detected into 
segments and characterizing the orientation of each segment. Each edge segment was 
then referenced to a map of expected proximal femur edge orientation and to a map of 
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the probability of edge location. Edge segments that did not conform to the expected 
orientation or which were in low probability regions were removed. Morphology 
operations were applied to the edge image(s) to connect any discontinuities. The edge 
image formed an enclosed boundary of the proximal femur. The region within the 
boundary was then combined with the binary image from global thresholding to form 
the final mask of the proximal femur. 

[0116] Automatic Segmentation and Measurement of the Femoral Cortex : Within 
a region of interest (ROI), edge detection was applied. Morphology operations were 
applied to connect edge discontinuities. Segments were formed within enclosed 
edges. The area and the major axis length of each segment were then measured. The 
regions were also superimposed on the original gray level image and average gray 
level within each region was measured. The cortex was identified as those segments 
connected to the boundary of the proximal femur mask with the greatest area, longest 
major axis length and a mean gray level about the average gray level of all enclosed 
segments within the proximal femur mask. 

[01 1 7] The segment identified as cortex was then skeletonized. The orientation of 
the cortex skeleton was verified to conform to the orientation map of the proximal 
femur edge. Euclidean distance transform was applied to the binary image of the 
segment. The values of distance transform value along the skeleton were sampled and 
their average, standard deviation, minimum, maximum and mod determined. 

[0118] Watershed Segmentation for Characterizing Trabecular Structure : 
Marrow spacing was characterized by determining watershed segmentation of gray 
level trabecular structures on the hip images, essentially as described in Russ "The 
Image Processing Handbook," 3 rd . ed. pp.494-501. This analysis takes the gray level 
contrast between the marrow spacing and adjacent trabecular structures into account. 
The segments of marrow spacing generated using watershed segmentation were 
measured for the area, eccentricity, orientation, and the average gray level on the x- 
ray image within the segment. Mean, standard deviation, minimum, maximum and 
mode were determined for each segment. In addition, various micro-structural and/or 
macro-anatomical parameters were assessed for several ROIs to predict the fracture 
path, as shown in Fig. 1 1. 
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[0119] Measurement of Femoral Neck BMD : DXA analysis of bone mineral 
density was performed in the femoral neck region of the femurs. 



[0120] Biomechanical Testing of Intact Femur : Each cadaveric femur sample 
(n=15) was tested for fracture load as follows. First, the femur was placed at a 15° 
angle of tilt and an 8° external rotation in an Instron 1331 Instrument (Instron, Inc.) 
and a load vector at the femoral head simulating single-leg stance was generated, 
essentially as described in Cheal et al. (1992) X Orthop. Res, 10(3):405-422. Second, 
varus/valgus and torsional resistive movements simulating passive knee ligament 
restraints were applied. Next, forces and movement at failure were measured using a 
six-degree of freedom load cell. Subsequently, a single ramp, axial compressive load 
was applied to the femoral head of each sample at 100 mm/s until fracture. Fracture 
load and resultant equilibrium forces and moments at the distal end of the femur were 
measured continuously. 

[0121] There was a weak positive correlation of femoral neck BMD (r=0.34, 
p=0.10) (Fig. 4A) and total BMD with failure load (r=0.28, p=0.15). Radiographs 
were analyzed in several regions of interest (ROI) at the femoral head, neck and 
proximal shaft to yield indices of trabecular micro-structure and macro-anatomic 
indices such as cortical thickness. The micro-structural parameter of Trabecular 
Segment Thickness from ROI 4 had the strongest failure load correlation coefficient, 
with r—0.75. Macro-anatomic indices such as Maximum Cortical Thickness of ROI 6 
and Median Cortical Thickness of ROI 5 correlate with failure load with r=0.65 
(p=0.005) and r=0.53 (p=0.02), respectively. 

[01 22] Based on these results, Trabecular Segment Thickness and Trabecular 
Separation from ROI 4 were combined to predict a failure load. Based on results 
from these 15 femora, correlation between predicted and actual failure loads was 
r=0.8 (pO.001) (Fig. 5B). The mean fracture load was 5.4 kiloNewton with a 
standard deviation of 2.3 kiloNewton. These statistics and the coefficients of 
multivariate linear regression were stored as data of the fracture load reference 
database. 
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[01 23] Influence of Positioning : The effects of the femur position were also 
examined in order to determine a set of measurements that are the most robust against 
the positioning variations that can occur during imaging. 

[01 24] Radiographs were taken at -1 5° (external rotation), and at every 5° 
increment up to +20° of internal rotation (70 kVp, photo-timer, centered on the 
femoral neck). Variability of a parameter was expressed as the coefficient of variation 
(COV) of the measurements at each angle. Of all the regions, ROI 4 showed the 
lowest average (root mean square) COV of the parameters. 

[01 25] As shown in Table 3 below, variability was generally lower for the 5° to 
15° range. This was also observed for the other regions of interest. Thus, internal 
rotation of 5° to 15° provides an acceptable margin of variation for a number of 
parameters. The regular AP hip x-ray imaging protocol used by technicians is 
therefore sufficient to control positioning variability. 
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Table 3 



Parameter 


Range of Rotation angles (degrees) 


0 10 

\J . . . M.XJ 


5 1 5 




1 o 9n 

lU...^iU 


1 ^ on 


Trab. Perimeter 

/TP a. 1 A 

/Total Area 


1 6 


0 s 

I/. J 


1 7 


1 0 




Free-end 
Inickness 


2.6 


0.5 


2.7 


3.3 


3.7 


Segment 
Tortuosity 


0.8 


1.0 


0.9 


0.7 


0.8 


Trabecular 

Segment 

Thickness 


2.5 


1.0 


3.6 


4.1 


3.0 


Trabecular 
Area Ratio 


2.5 


1.2 


6.9 


7.7 


5.3 


Trabecular 
Bone Pattern 
Factor 


1.3 


1.3 


4.1 


4.4 


3.8 


Trabecular 
Separation 


1.9 


2.0 


1.6 


0.6 


1.8 


Segment 
Solidity 


4.1 


2.6 


4.7 


5.4 


3.9 



[0126] Influence of Radiographic Exposure Settings : The influence of hip x-ray 
exposure variations on image quality and the subsequent analysis of structural 
parameters were also tested. 

[01 27] The right hip of a frozen cadaver pelvis was imaged with 60kVp, 70kVp, 
and SOkVp at 150mA with automatic exposure using the photo-timer, followed by an 
exposure one step below and another at one step above the auto exposure, in terms of 
niAs. An additional image was taken at 75kVp with the photo timer. 

[01 28] Most parameters exhibited the least variation (across mAs) at 60kVp ? with 
variability growing with increased kVp. Trabecular separation measurements in ROI 
7 had COV's of 2.1% , 4.2% and 5.1% at 60kVp, 70kVp and 80kVp, respectively. 
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These may represent the variability to be expected when using manual time settings in 
the absence of the automatic phototimer function. 



[01 29] When measurements from only the images captured using the automatic 
phototimer were considered across kVp (60, 70, 75, 80), the most reproducible 
measurement was trabecular perimeter, with an average GOV of 1.9%. 

[01 30] Our results indicate -that the use of a phototimer can markedly reduce the 
variability of exposures due to subjective kVp setting choices. Radiographs produced 
with proper and consistent use of the phototimer had acceptable variations of micro- 
structural and macro-anatomical measurements. 

Example 2: Correlation of 2D and 3D measurements 

[0131] To demonstrate that these methods that use 2D x-ray technology to 
quantitatively assess trabecular architecture are as effective as 3D /aCT, which serves 
as a gold standard for such measurements, the following experiments were performed. 
Bone cores (n=48) were harvested from cadaveric proximal femora. Specimen 
radiographs were obtained and 2D structural parameters were measured on the 
radiographs. Cores were then subjected to 3D p,CT and biomechanical testing. The 
\xCT images were analyzed to obtained 3D micro-structural measurements. Digitized 
2D x-ray images of these cores were also analyzed as described herein to obtain 
comparative micro-structural measurements. 

[01 32] Results showed very good correlation among the numerous 2D parameters 
and 3D juCT measurements, including for example correlation between 2D Trabecular 
Perimeter/Trabecular Area (Tb.P/Tb.A) with 3D Bone Surface/Bone Volume (r = 
0.92, p< 0.001), and 2D Trabecular Separation (Tb.Sp) with 3D Trabecular Separation 
(r = 0.88, p < 0.001), as shown in Fig. 3. The 2D Tb.P/Tb.A and 2D Tb.Sp also 
function correlate very well as predictive parameters for the mechanical loads 
required to fracture the cores, with r = -0.84 (p <0.001) and r = -0.83 (p < 0.001), 
respectively, when logarithmic and exponential transformations were used in the 
regression, as shown in Fig. 4. 
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[0133] These results demonstrate that 2D micro-structural measurements of 
trabecular bone from digitized radiographs are highly correlated with 3D 
measurements obtained from \iCT images. Therefore, the mechanical characteristics 
of trabecular bone micro-structure from digitized radiographic images can be 
accurately determined from 2D images. 

Example 3: Sliding Window Analysis and Watershed 

Segmentation 

[01 34] To show feasibility of the approach to better predict an individual's failure 
load and fracture risk, a sliding window analysis of the proximal femur in 3 cadaveric 
samples was also performed. Instead of using fixed ROI's as described in Fig. 2, a 
regular grid was laid over the proximal femur in the x-ray taken before the mechanical 
failure tests. A rectangular region of interest of a fixed size was placed at each grid 
node, and bone structure parameters were evaluated within the boundaries of the ROI 
at each position. This resulted in a value for each bone parameter at each grid node, 
which was displayed in a color-coded map of the proximal femur for each parameter. 

[0135] The pre-fracture x-rays were then aligned with, the post-fracture images, so 
that the fracture lines can be shown in the color maps, as shown in Fig. 6. It can be 
seen in the samples presented in Fig. 6 that certain parameters (e.g., for Area Ratio 
and Trabecular Perimeter ) have a very good agreement between low values 
("valleys") in the color maps and the fracture lines, suggesting that a sliding window 
bone structural analysis can be used to generate a prediction of the exact location 
where the bone will fracture. The valleys can be found by applying a watershed 
transformation to the negative values of the parameter maps. For other bone 
parameters that indicate bone weakness with high values, the watershed 
transformation can be applied directly to the map. 

Example 4: Prediction of fracture risk using fracture load 
reference database. 

[01 36] A hip x-ray of a cadaver pelvis was exposed using standard clinical 
procedure and equipment. The radiograph film was developed and digitized. The 
image was then analyzed to obtain micro-structure, and macro-anatomical parameters. 
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The local maximum spacing, standard deviation of cortical thickness of ROI 3, 
maximum cortical thickness of ROI 5, and mean node-free end length for ROI 3 were 
used to predict load required to fracture the cadaver hip using the coefficients of 
multivariate linear regression stored in the fracture load reference database. The 
predicted fracture load was 7.5 kiloNewton. This fracture load is 0,98 standard 
deviation above the average of the fracture load reference database (or z-score = 
0.98). This result may suggest that the subject had a relatively low risk of sustaining a 
hip fracture as compared to the population of the reference database. 

Example 5. Prediction of hip fracture risk from radiographic 

Images 

[0137] Individualized hip fracture risk is determined as shown in Fig. 7 A. 
Briefly, an x-ray of the hip is taken at step 70 L 

[0138] At step 702, a 2D fracture line is predicted. The micro- and macro- 
architecture of the proximal femur in the image is determined by performing 
automated analyses, as described in Examples 1 and 4. Algorithms for analysis of 
density, length, thickness, and orientation of trabeculae as well as cortical bone 
thickness in an ROI in the radiograph are developed using Matlab (The MathWorks, 
Inc., Natick, MA). Similarly, software is developed for automated sliding window 
analysis (Example 3) of parameters, including at least one bone structure parameter, 
to produce a distribution map of the proximal femur for each parameter. 

[01 39] In addition, local abnormalities of bone structure are determined from the 
parameters maps generated as described in Example 3. Regions of high or low values 
will be evaluated to determine bone strength patterns and used to predict a location of 
hip fracture. 

[0140] The parameter maps generated as described in Example 3 are used to 
identify regions on the bone that have abnormal local structural properties. These 
regions of high or low values can indicate patterns for stronger or weaker 
characteristics of bone. The parameter maps generated from the hip x-ray provide a 
spectrum of trabecular characteristics and can be interpreted as spatial distributions of 
bone strength. They will be used to predict the location of a hip fracture. 
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[0141] In a first step, depending on the kind of parameter, the low values 
("valleys") or high values ("ridges") on the parameter maps are traced using a 
watershed segmentation operation (see Fig. 8), The resulting boundaries between the 
regions are regarded as potential fracture lines. 

[0142] In a second step, the path along the potential fracture lines that is most 
likely to coincide with the actual fracture location is determined. A two-dimensional 
curved beam model as described by Mourtada et al. (1996) J Orthop Res, 14(3):483- 
492 is applied on the thresholded x-ray image of the proximal femur. First, the 
femoral shaft and neck are approximated by linear axes. The axis in the 
intertrochanteric region is approximated by a hyperbolic curve that is asymptotic to 
the linear axes of the femoral neck and shaft, with the focus point at the neck-shaft 
angle bisector (Fig. 9). Given a loading condition, the internal bending moment, M 9 is 
calculated along the neutral axis at 1 mm intervals. For regions where curvature is 
negligible, the normal stresses along the boundary of the femur are calculated by 
Equation (1): 

M'X n . 

<T n = (1) 

n CSMI 

where CSMI is cross-sectional moment of inertia, and x is the distance from the 
neutral axis to the point where stress is calculated. Since peak stresses occur on the 
surfaces, x will be the perpendicular distance of the surface boundary to the neutral 
axis. Along the curvature, normal stresses can be calculated using Equation (2), 

M • x 

cr = (2) 

n CSAe-^-x) 

where CSA is the cross-sectional area, e is the distance between the centroid axis and 
neutral axis, and R na is the radius of curvature of the neutral axis. The loading 
condition applied to the curved beam model will simulate a fall on the side from 
standing height with the estimated forces obtained using the methods below. Both 
CSMI and CSA can be estimated for each cross section by integrating the optical 
density over the section profile. The relative density and stress values are sufficient 
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for the purpose of locating fracture location, 
be insignificant over the proximal femur. 
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[0143] Two common types of fracture, intertrochanteric and femoral neck, will be 
considered. Two tensile peaks are known to exist on the medial surface for the fall 
loading condition (see, e.g., Mourtada et aL, supra). The peak closer to the bisector 
of the neck-shaft angle is identified as the starting point of the intertrochanteric 
fracture, and the other one that is known to be on the posterior surface of the neck, as 
the starting point for the femoral neck fracture. The cross-sectional lines 
corresponding to the position of tensile peaks will be considered. The predicted 
fracture paths will be traced by selecting the watershed boundary segments that are 
closest to the selected cross-sectional lines (Fig. 10). 

[01 44] To predict the likelihood of intertrochanteric or femoral neck fracture, the 
values of the parameter map underlying the selected fracture paths will be evaluated 
and compared. The more likely fracture path will have a lower mean value of a 
parameter that represents bone strength as optimized by cadaver mechanical loading 
tests. 

[0145] At step 703, a local micro- and macro-structural analysis along the 
predicted fracture line is performed to estimate the load at which the bone will 
fracture in a particular falling scenario. The case-specific ROI is placed around the 
predicted fracture line in the trabecular bone using a region growing technique with 
value constraints. Cortical bone parameters are evaluated in the areas adjacent to the 
trabecular bone ROI with boundaries determined by perpendicular projection of the 
outer contact points between trabecular ROI and cortical bone onto the outer contour 
of the cortical bone (Fig. 11). Multivariate regression will then be used to calculate a 
failure load Ff a ii U re from the results of the different bone parameter analyses. 

[01 46] The risk of sustaining an osteoporotic hip fracture does not only depend on 
the femoral failure load, but also on the impact on the femur in a fall. Factors that 
influence the severity of the impact are, among others, soft tissue thickness, standing 
height, and body mass. Impact decreases with increasing soft tissue covering the hip, 
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while it is increased with greater standing height or body mass, 
calculated at step 704. 
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[0147] Estimation of femoral impact is performed essentially as described in 
Kroonenberg et al. (1995) J. Biomech Eng. 117(3):309-318. Calculations shown in 
Equations. (3) - (8) below are based on studies of women. Equations for men can be 
derived accordingly. 

[01 48] The hip impact velocity is given by Equation (3) 5 
V = 2.72Vh (3), 

where h is the full body height. The effective mass, i.e., the mass of the part of the 
body that contributed to the impact force on the hip is calculated as shown in 
Equation (4): 

7 

M = — -m (4) 
where m is the total body mass. 

[01 49] The peak force on the greater trochanter can then be approximated by 
Equation (5), at step 705: 

= vV^ = L6Vh^k (5) 

[01 50] The soft tissue stiffness k correlates negatively with soft tissue thickness 
(see, Robinovitch et al. (1991) J. Biomech Eng. 1 13:366-374). Fitting the data 
obtained by Robinovitch et al. for loads with 100 N with a power curve, Equation (6) 
is obtained: 

k = 486x" 083 (6) 

with x being the soft tissue thickness. Soft tissue stiffness dependency on loading 
force can be approximated by Equation (7): 

-F 

k = 71060- (1-e 151 ) (7) 
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[0151] Using Equation. (7) for 100 N, k= 34415 N/m for women. Since with 
Equation (7) soft tissue stiffness plateaus at 71 kN/m for loads expected in a fall, 
Equation (6) by 71000 / 34415 = 2.1 to obtain the relationship between soft tissue 
thickness and soft tissue stiffness for these higher loads as given by Equation (8): 

k = 102hT 0 ' 83 (8) 

[0152] The soft tissue thickness x will be measured in the hip x-ray laterally 
between the greater trochanter and the skin. 

[0153] As discussed above with Example 4, a fracture load reference database can 
be used for more accurate determination of the fracture load. 

[01 54] At step 7 10, a measure of fracture risk can then be calculated as the ratio 
of the peak impact force obtained at step 705 via Equation (5) and the predicted 
failure load obtained at step 703: 

IFRI (Individualized Fracture Risk Index) = F pe ak/Ff a ii urc (9) 

[01 55] Thus, when the Individualized Fracture Risk Index is low (IFRI «1), the 
forces applied to the femur are much lower than required to fracture it, and the bone is 
at low risk of failure. However, when the IFRI is high (IFRI »1), failure of the bone 
is predicted to occur. 

[01 56] Fig. 7B summarizes the procedure for predicting hip failure load from an 
x-ray of the proximal femur. As shown, step 702 for predicting fracture lines includes 
two sub-steps 7021 and 7022. Sub-step 7021 uses sliding ROI analysis to analyze 
bone parameters, including at least one bone structure parameter, for each sampling 
point and generate bone parameter maps. Sub-step 7022 uses watershed segmentation 
of parameter maps to identify most likely fracture lines. Step 703 for predicting the 
failure load includes sub-steps 7031 and 7032. Sub-step 7031 places ROI along 
predicted fracture line and analyze trabecular and cortical bone parameters; and sub- 
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step 7032 combines multiple bone parameters and calculates fracture load using 
multivariate regression. 



[0157] Fig. 12 is a block diagram of a computer program used to predict fracture 
risk according to one embodiment of the present invention. As shown, a module 1201 
receives images of skeleton parts of patients and/or references. A module 1202 
receives the images from the module 1201 and derives bone structure parameters 
therefrom. A module 1203 measures a fracture load of a skeleton part of a reference. 
A module 1204 correlates the reference's fracture load to the reference's bone 
structure parameter, taking clinical risk factors of the reference into consideration. A 
module 1205 controls storage of the correlation. A module 1206 generates a 
parameter map from the derived bone structure parameter of a patient to predict a 
fracture line. A module 1207 receives the correlation from the module 1204 and the 
fracture line from the module 1206, and calculates the fracture load of the patient. A 
module 1208 estimates a body habitus of the patient. A module 1209 receives the 
body habitus estimation from the module 1208 and the fracture load of the patient 
from module 1207, and calculates a peak input force. The risk of fracture is 
predicted at a module 1210 by calculating the ratio between the fracture load of the 
patient from the module 1207 and the peak impact force from the module 1209. 

[01 58] Although Fig. 12 illustrates modules of a computer program, a skilled 
artisan would appreciate that hardware and firmware could be used to realize the 
functions of the modules, and the modules could be distributed at different locations. 
A suitably programmed computer can constitute hardware counterparts of each of the 
modules in Fig. 12. 

[01 59] The foregoing description of embodiments of the present invention has 
been provided for the purposes of illustration and description. It is not intended to be 
exhaustive or to limit the invention to the precise forms disclosed. Many 
modifications and variations will be apparent to the practitioner skilled in the art. The 
embodiments were chosen and described in order to best explain the principles of the 
invention and its practical application, thereby enabling others skilled in the art to 
understand the invention and the various embodiments and with various modifications 
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that are suited to the particular use contemplated. It is intended that the scope of the 
invention be defined by the following claims and their equivalents. 
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CLAIMS 

1. A method for diagnosing a musculoskeletal disease, comprising: 
obtaining an image of a part of a skeleton of a target; 
locating at least one region of interest on the image of the target; 
extracting image data from the image of the target; 

deriving at least one bone structure parameter from the image data of the 
.target; and 

predicting a fracture with the bone structure parameter of the target. 

2. The method of claim 1, further comprising: 

obtaining an image of the same part of the skeleton of a reference; 
locating at least one region of interest on the image of the reference; 
extracting image data from the image of the reference; 
deriving at least one bone structure parameter from the image data of the 
reference; 

measuring a fracture load of the skeleton part of the reference; and 
correlating the reference's bone structure parameter to the reference's fracture 

load. 

3. The method of claim 2, further comprising predicting the fracture load of the 
target with the target's bone structure parameter and the correlation of the reference's 
bone structure parameter and the reference's fracture load. 

4. The method of claim 2, further comprising: 

creating a fracture load reference database with data about at least one bone 
structure parameter and fracture load of skeleton parts of at least two references; and 
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predicting the fracture load of the target with the target's bone structure 
parameter and data from the fracture load reference database. 

5. The method of claim 4, wherein the data is stored based on clinical risk factors 
of the references. 

6. The method of claim 4, further comprising creating a fracture risk model 
correlated with fracture load by statistically analyzing at least one bone structure 
parameter of the references. 

7. The method of claim 6 5 wherein the fracture risk model is a multivariate linear 
model derived from multivariate linear regression of at least one bone structure 
parameter and clinical risk factors of the references. 

8. The method of claim % further comprising deriving a second parameter from 
the image data of the target and the reference. 

9. The method of claim 8, wherein the second parameter is related to bone 
mineral density. 

10. The method of claim 8, wherein the second parameter is related to bone 
biomechanical characteristics. 

1 1 . The method of claim 8 5 wherein the second parameter is related to bone 
structure. 

12. The method of claim 2, further comprising generating a parameter map from 
the derived bone structure parameter of the target to predict a fracture line. 

13. The method of claim 12, wherein the parameter map is derived using 
statistical comparisons of the derived bone structure parameter to a reference 
population. 

14. The method of claim 12, wherein the parameter map is generated via a sliding 
window region of interest analysis of the derived bone structure parameter. 
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15. The method of claim 12, further comprising identifying local abnormalities of 
bone structure from the parameter map. 

16. The method of claim 12, further comprising: 

tracing low values or high values on the parameter map; and 
determining a potential fracture line from the low values or high values. 

17. The method of claim 12, further comprising using watershed segmentation of 
parameter maps to identify the fracture line, 

18. The method of claim 1, further comprising generating a parameter map from 
the derived bone structure parameter of the target to predict a fracture line. 

19. The method of claim 12, further comprising analyzing the at least one bone 
structure parameter along the predicted fracture line to predict a fracture load at which 
a fracture will occur. j 

20. The method of claim 12, further comprising: 

deriving a second bone parameter from the image of the target and the 
reference; 

combining the bone structure parameter and the second bone parameter; and 
calculating the fracture load using multivariate regression. 

21 . The method of claim 1 8, further comprising analyzing the at least one bone 
structure parameter along the predicted fracture line to predict a load at which a 
fracture will happen. 

22. The method according to claim 19, further comprising estimating a body 
habitus of the target, 

23. The method according to claim 22, wherein the body habitus is related to a 
soft tissue thickness of the target. 
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24. The method according to claim 22, wherein the body habitus is related to a 
standing height of the target. 

25. The method according to claim 22, wherein the body habitus is related to a 
body mass of the target. 

26. The method according to claim 1, further comprising estimating the body 
habitus of the target. 

27. The method of claim 22, further comprising calculating a peak impact force on 
the skeleton part when the target falls. 

28. The method of claim 27, further comprising predicting a fracture risk by 
calculating the ratio between the peak impact force and the fracture load. 

29. The method of claim 1, wherein the bone structure parameter is a bone micro- 
structure parameter. 

30. The method of claim 1, wherein the bone structure parameter is a bone macro- 
structure parameter. 

3 1 . The method of claim 1, wherein the at least one region of interest is located 
automatically. 

32. The method of claim 1, wherein the image is selected from the group 
consisting of x-ray images, x-ray tomosynthesis, ultrasound, computed tomography, 
magnetic resonance imaging, optical coherence tomography, single photon emission 
tomography, and positron emission tomography. 

33. The method of claim 1, wherein the image is a 2D image. 

34. The method of claim 1, wherein the image is a 3D image. 

35. The method of claim 1, wherein the image is a 4D image. 

36. The method of claim 1, further comprising converting the image to an image 
of more dimensions. 
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37. The method of claim 1, further comprising repeating the step for obtaining an 
image. 

38. The method of claim 1, further comprising repeating the step for locating at 
least one region of interest. 

39. The method of claim 1, further comprising repeating the step for extracting 
image data from the image. 

40. The method of claim 1, further comprising repeating the step for deriving at 
least one bone structure parameter. 

41 . The method of claim 1, further comprising: 
transmitting the image to a second location; 

converting the image to a pattern of normal or diseased using the bone 
structure parameter; and 

analyzing the converted image. 

42. The method of claim 41, further comprising transmitting the pattern to a third 
location for analyzing. 

43. A database storing: 

data about at least one bone structure parameter and fracture load of skeleton 
parts of at least two references obtained according to the method of claim 19. 

44. The database of claim 43, wherein the data is stored based on clinical risk 
factors of the references. 

45. The database of claim 44, wherein the clinical risk factor is the reference's 
age. 

46. The database of claim 44, wherein the clinical risk factor is the reference's 
gender. 
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47. The database of claim 44, wherein the clinical risk factor is the reference's 
weight. 

48. The database of claim 43 , further storing a fracture risk model correlated with 
fracture load. 

49. The database of claim 44, further storing a fracture risk model correlated with 
fracture load. 

50. The database of claim 48 ? wherein the fracture risk model is generated by 
statistically analyzing at least one bone structure parameter of the references. 

5 1 . The database of claim 49, wherein the fracture risk model is generated by 
statistically analyzing at least one bone structure parameter of the references. 

52. The database of claim 50, wherein the fracture risk model is coefficients of a 
multivariate linear model derived from multivariate linear regression of at least one 
bone structure parameter and clinical risk factors of the references. 

53. A method for monitoring a musculoskeletal disease, comprising: 
diagnosing a fracture risk of a target at a first time according to the method of 

claim 1; 

diagnosing the fracture risk of a target at a second time according to the 
method of claim 1 ; 

comparing the target's fracture risk at the first time and at the second time. 

54. A method for determining effect of a treatment on a target's musculoskeletal 
disease, comprising: 

diagnosing a fracture risk of the target at a first time according to the method 
of claim 1; 

providing the treatment to the target; 
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diagnosing the fracture risk of the target at a second time according to the 
method of claim 1; and 

comparing the target's fracture risk at the first time and the second time. 

55. The method of claim 54, wherein the treatment is administering a candidate 
agent. 

56. A computer program product containing program code for performing the 
method according to claim 1 . 

57. A computer program product containing program code for performing the 
method according to claim 2. ; 

58. A computer program product containing program code for performing the 
method according to claim 3. 

59. A computer program product containing program code for performing the 
method according to claim 4. 

60. A computer program product containing program code for performing the 
method according to claim 5. 

61 . A computer program product containing program code for performing the 
method according to claim 6. 

62. A computer program product containing program code for performing the 
method according to claim 12. 

63. A computer program product containing program code for performing the 
method according to claim 14. 

64. A computer program product containing program code for performing the 
method according to claim 17. 

65. A computer program product containing program code for performing the 
method according to claim 19. 
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66. A computer program product containing program code for performing the 
method according to claim 20. 

67. A computer program product containing program code for performing the 
method according to claim 22. 

68. A computer program product containing program code for performing the 
method according to claim 27. 

69. A system for analyzing musculoskeletal-related data, comprising: 
means for receiving an image of a part of a skeleton; 

means for deriving at least one bone structure parameter from the image; and 
means for calculating a possibility of a fracture of a target with a bone 
structure parameter of the target. 

70. The system according to claim 69, further comprising: 

means for obtaining a fracture load of a skeleton part of a reference; and 
means for correlating a bone structure parameter of the reference to the 
fracture load of the reference. 

71 . The system according to claim 70, further comprising means for receiving the 
target's bone structure parameter and the correlation of the reference's bone structure 
parameter and the reference's fracture load, and calculating the fracture load of the 
target. 

72. The system according to claim 7 1 further comprising means for storing the 
correlation of the reference's bone structure parameter and the reference's fracture 
load. 

73. The system according to claim 72, wherein the means for storing also receives 
clinical risk factors of the reference. 
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74. The system according to claim 71, further comprising means for generating a 
parameter map from the target's bone structure parameter for presenting a possible 
fracture line. 

75. The system according to claim 74, further comprising means for analyzing the 
target's bone structure parameter along the possible fracture line to calculate a 
fracture load at which a fracture will occur. 

76. The system of claim 75 5 further comprising means for estimating a body 
habitus of the target. 

77. The system of claim 76, wherein the body habitus is related to a soft tissue 
thickness of the target. 

78. The system of claim 76, wherein the body habitus is related to a standing 
height of the target. 

79. The system of claim 76, wherein the body habitus is related to a body mass of 
the target. 

80. The system of claim 76, further comprising means for calculating a peak 
impact force on the skeleton part when the target falls. 

8 1 . The system of claim 80, further comprising means for calculating the ratio 
between the peak impact force and the fracture load of the target. 

82. The system of claim 69, wherein the bone structure parameter is a bone 
microstructure parameter. 

83. The system of claim 69, wherein the bone structure parameter is a bone 
macro-structure parameter. 
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